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ABSTRACT 

We measure a tilt of 86 ± 6° between the sky projections of the rotation axis of the WASP-7 star, and the 
orbital axis of its close-in giant planet. This measurement is based on observations of the Rossiter-McLaughlin 
(RM) effect with the Planet Finder Spectrograph on the Magellan II telescope. The result conforms with the 
previously noted pattern among hot-Jupiter hosts, namely, that the hosts lacking thick convective envelopes 
have high obliquities. Because the planet's trajectory crosses a wide range of stellar latitudes, observations 
of the RM effect can in principle reveal the stellar differential rotation profile; however, with the present data 
the signal of differential rotation could not be detected. The host star is found to exhibit radial-velocity noise 
("stellar jitter") with an amplitude of w30 m s"' over a timescale of days. 

Subject headings: techniques: spectroscopic — stars: rotation — planetary systems — planets and satellites: 
formation — planet-star interactions — stars: individual (WASP-7) 



1. INTRODUCTION 

In the solar system the Sun's equatorial plane is aligned to 
within 7° with the ecliptic. For many stars in systems with 
close-in gas giants ("hot-Jupiters"), this is not the case. Over 
the last three years it was found that some hot Jupiters have 
highly inclined or even retrograde orbits with respect to the ro- 
tational spins of their host stars (see e.g. Hebrard et al. 2008; 
Winn et al. 2009; Johnson et al. 2009; Triaud et al. 2010; Winn 
et al. 2011; Simpson et al. 2011). Understanding what causes 
these orbital tilts and why some hot Jupiters are well-aligned 
with their parent stars might aid our understanding of why 
these giant planets are found so close to their host stars, com- 
pared to Jupiter. 

Different classes of processes have been proposed which 
might transport giant planets from their presumed birthplaces 
at distances of many astronomical units from their host stars, 
inward to a fraction of an astronomical unit, where we find 
them. Some of these processes are expected to change the 
relative orientation between the stellar and orbital spin (e.g. 
Nagasawa et al. 2008; Fabry cky & Tremaine 2007; Chatterjee 
et al. 2008), while others will conserve the relative orientation 
between orbital and stellar spin (e.g. Lin et al. 1996), or even 
reduce a misalignment between them (Cresswell et al. 2007). 

Winn et al. (2010) and Schlaufman (2010) found that close 
in giant planets tend to have orbits aligned with the stellar 
spin if the effective temperature (Feff) of their host star is 
< 6250 K and misaligned otherwise. Winn et al. (2010) fur- 
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ther speculated that this might indicate that all giant planets 
are transported inward by processes that create large obliqui- 
ties. In this picture, tidal torques exerted on the star by the 
close in planet realign the two angular momentum vectors. 
The realignment time scale would be short for planets around 
stars with convective envelopes (Teff < 6250K), but long if the 
star does not have a convective envelope (> 6250K). Adding 
to this picture, Triaud (2011) recently argued that relatively 
young stars have high obliquities, while older stars are ob- 
served to have low obliquities. 

Here we present measurements of the spin-orbit angle in 
the WASP-7 system. The planet WASP-7b was discovered by 
Hellier et al. (2009) and found to have a mass of 0.96 Mjup. 
The host star has a mass of 1.28M0, a projected rotation 
speed of 17 ±2 kms"' (Hellier et al. 2009), an effective tem- 
perature of 6520 ± 70 K, and a solar metallicity ([Fe/H] = 
0.0 ±0.1) (Maxted et al. 201 1). Based on the aforementioned 
pattern, we would expect that our measurement would show a 
misalignment between orbital and stellar spins. 

This article is organized as follows. The following section 
describes the new spectroscopic data, and the analysis of the 
Rossiter-McLaughlin effect. Section 4 considers some possi- 
ble explanations for the high level of noise in the radial veloc- 
ities, including the possibility of an eccentric orbit. Section 5 
discusses the impact of the radial-velocity noise on the mea- 
surement of the stellar obliquity. This section also presents an 
attempt to detect the differential rotation of the host star using 
the Rossiter-McLaughlin effect. 

2. OBSERVATIONS 

We observed WASP-7 with the Magellan II (Clay) 6.5 m 
telescope and the Planet Finder Spectrograph (PES; Crane 
et al. 2010). We gathered 37 spectra spanning the transit of 
2010 August 27/28. The integration times were lOmin and 
the complete sequence spanned ^7.5 hr The stellar spec- 
tra were observed through an iodine gas cell, imprinting a 
dense forest of sharp absorption lines on the stellar spectra to 
help establish the wavelength scale and instrumental profile. 
During the transit night, an additional spectrum was obtained 
without the iodine cell, to serve as a template spectrum for 



TABLE 1 

Relative Radial Velocity measurements of WASP-7 




Fig. 1. — Radial velocities of WASP-7 before, during, and after tlie 
transit of its planet. The radial velocities are plotted as a function of time 
from inferior conjunction. The upper panel shows the measured RVs and the 
best-fitting model. The dashed line shows the same RM model, but with an 
orbital model with parameters fixed at those presented by Hellier et al. (2009) 
(see also section 4). In the middle panel, the apparent orbital contribution to 
the observed RVs has been subtracted, thereby isolating the RM effect. The 
lower panel shows the residuals. The light and dark gray bars in the lowest 
panel indicate times of first, second, third, and fourth contact. 

relative radial-velocity (RV) determination. 

3. ANALYSIS OF THE ROSSITER-MCLAUGHLIN EFFECT 

To derive the relative RVs we compared the spectra ob- 
served through the iodine cell with the stellar template spec- 
trum multiplied by an iodine template spectrum. The velocity 
shift of the stellar template as well as the parameters of the 
point-spread function (PSF) of the spectrograph are free pa- 
rameters in this comparison. The velocity shift of the template 
that gives the best fit to an observed spectrum represents the 
measured relative RV. In particular we used a code based on 
that of Butler et al. (1996). The RVs are presented in Table 1 
and displayed in Figure 1 . 

We take advantage of the Rossiter-McLaughlin (RM) effect 
to measure the projected angle between the orbital and stel- 
lar spins (A) and the projected stellar rotation speed (vsin/*). 
Here v indicates the stellar rotation speed and i^, the inclina- 
tion of the stellar spin axis towards the observer. The RM ef- 
fect is a spectroscopic distortion of the rotationally-broadened 
stellar absorption lines, which occurs when a companion star 
or planet is in front of the star and hides part of the rotating 
stellar surface from the observer's view. The position of the 
distoition on the stellar absorption line depends on the radial 
velocity of the hidden portion of the stellar surface. Therefore 
the observed distortion of the stellar absorption lines can be 
connected to the geometry of the transit. This shape change 
can be measured directly and relevant parameters can be de- 
rived (Albrecht et al. 2007; Collier Cameron et al. 2010). 

Our analysis is divided into 7 parts. In section 3.1 we dis- 
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cuss the data qualitatively. Sections 3.2-3.6 discuss different 
phenomena which affect the shape of stellar absorption lines. 
Section 3.7 discusses the model of the RM effect which we 
adopted, and presents the quantitative results. 

3.1. Qualitative expectations 

The RM effect is evident in Figure 1 as the large negative 
velocity excursion (blueshift) that was observed throughout 
the transit. The effect was observed with a high signal-to- 
noise ratio. Simply from the observation that the effect is a 
blueshift throughout the transit, we may obtain some informa- 
tion on the spin-orbit alignment of the system. A qualitative 
discussion will help in understanding the quantitative analysis 
to be discussed later in this paper. 

If the projections of the stellar and orbital spins were 
aligned, then the planet would first traverse the half of the 
stellar surface for which the rotation velocity has a component 
directed toward the observer (blueshifted). The blockage of a 
portion of this blueshifted half of the star would cause the ab- 
sorption lines to appear slightly redshifted. Then, in the sec- 
ond half of the transit, the reverse would be true: the anoma- 



The rotation of the oblique star WASP-7 



3 




-30 -20 -10 10 20 30 -2 -1 1 2 

velocity [km s"^] time [lir] 



Fig. 2. — Line broadening mechanisms and their influence on the RM effect. Left. — Model absorption line exhibiting the RM effect. The blue dashed line 
shows the idealized case in which the only line broadening mechanism is uniform rotation. The red dash-dotted line shows the case in which macroturbulence 
and the convective blueshift are also taken into account. The black line shows the case in which the former effects as well as differential rotation are taken into 
account. The lower set of lines indicate the "loss of light" in the stellar absorption line due to shadowing by the planet. For visibility this was multiplied by a 
factor of 3. These models are based on the parameters given in Table 2. Right. — Time variation of the RM effect (anomalous RV) for the same models as in the 
left panel. The circles indicate the particular phase for which the model line profiles in the left panel are shown. While the shift in the first moment appears to be 
larger for the blue dashed line (left panel) than for the other lines, its RM effect is smaller (right panel). This is partly a result of the measurement technique as 
described in section 3.5. In addition, the black and red dash-dotted lines are influenced by the convective blueshitt (section 3.3). 



lous RV would be a blueshift. In contrast, the RM effect was 
observed to be a blueshift throughout the transit. This implies 
that the planet's trajectory is entirely over the redshifted half 
of the star. 

The symmetry of the RM signal about the midtransit time 
provides further information. For a star in uniform rota- 
tion, the line-of-sight component of the rotation velocity at 
any point on the photosphere is proportional to the distance 
from the projected rotation axis; see, e.g., p. 461-462 of Gray 
(2005). Therefore if the planet's trajectory is at a constant 
distance from the projected rotation axis (i.e. if the projected 
orbital and spin axes are perpendicular), the hidden velocity 
component will be nearly constant in time and the only vari- 
ation in the RM signal arises from limb darkening, which is 
symmetric about the midtransit time. In fact this is the only 
way to produce a time-symmetric RM effect unless the impact 
parameter is nearly zero, which is not the case for WASP-7 
(Southworth et al. 201 1). 

We conclude that the projected orbital and spin axes are 
nearly perpendicular, and since the anomalous RV is a 
blueshift we expect A « 90° (as opposed to -90°), using the 
coordinate system of Ohta et al. (2005). This quantitative con- 
clusion is confirmed by our quantitative analysis discussed be- 
low; but first we must take into account several effects besides 
rotation that influence the shape of the measured absorption 
lines. 

3.2. Line broadening and the impact on the analysis 

Stellar absorption lines are broadened not only by the 
Doppler shifts due to stellar rotation, but also by random mo- 
tions of material on the visible stellar surface. These ran- 
dom motions are often referred to as microturbulence and 
macroturbulence, depending on whether the length scale of 
the velocity field is smaller or larger than the length scale 
over which the optical depth is of order unity (see, e.g. Gray 



2005). In this picture, the influence of microturbulence can 
be described by a simple convolution by a Gaussian function 
of an appropriate width with the rotationally-broadened stel- 
lar lines. The effect of the macroturbulence depends on the 
angle between the line of sight and the local normal to the 
stellar surface. Near the center of the stellar disk, the Doppler 
shifts are produced mainly by motions that are radial with re- 
spect to the stellar center; whereas near the stellar limb, the 
Doppler shifts are produced by motions tangential to the stel- 
lar surface. Therefore modeling this effect requires a spatial 
integration over the visible hemisphere. For this purpose we 
employ the semi-analytical approach of Hirano et al. (201 lb). 

Other contributions to the line width come from collisional 
broadening, and Zeeman splitting. These effects are small for 
WASP-7 compared to the other effects, and we describe them 
by the convolution of the disk-integrated line profiles with a 
Lorentzian function of width 1 km s"' . 

3.3. Convective blueshift 

Another potentially relevant effect on line profiles is the 
convective blueshift. At the top of a convective cell, where 
it meets the photosphere, the rising material is hotter and 
therefore more luminous than the sinking material. Due to 
the dominance of the rising material over the falling material, 
an observer placed vertically above the stellar surface would 
detect an overall net blueshift in the integrated light. This 
effectively outward radial velocity is referred to as the con- 
vective blueshift. When integrating over the visible stellar 
disk, the convective blueshift from the center of the disk has 
a maximal effect on the observed radial velocity, while the 
convective blueshift from points near the stellar limb have a 
weaker effect because the outward radial velocity has only a 
small component along the line of sight. Thus, the influence 
of the convective blueshift is strongest near the center of the 
line profile, and weaker in the wings of the lines where a sub- 
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Fig. 3. — Photometry of WASP-7 transits. The upper panel shows the 
light curve obtained by Southworth et al. (201 1) in the Gunn / filter, and our 
best-fitting model. The lower panel shows the residuals between the data and 
best-fitting model. 



stantial portion of the light originates from the stellar limb. 

Therefore, the disk-integrated light has not only an overall 
Doppler shift but also an asymmetry in the stellar absorption 
lines. To describe this, we employ the model of Shporer & 
Brown (2011). This model captures the first-order effects of 
the convective blueshift but ignores a higher-order asymme- 
try in the absorption lines: the line cores are formed relatively 
high in the photosphere where turbulent motions may be less 
important. As WASP-7 is a relatively hot star, with high turbu- 
lent velocities, it seems appropriate to ignore the latter effect 
(see Figure 17.15 in Gray 2005). 

3.4. Differential rotation 

When analyzing the RM effect researchers normally as- 
sume the eclipsed star is rotating uniformly (no differential 
rotation). In the case of WASP-7, the planet's trajectory is ev- 
idently perpendicular to the projected stellar equator (see sec- 
tion 3.1), and therefore spans a wide range of stellar latitudes, 
and the effects of differential rotation might be expected to be 
especially important. Ignoring the possibility of differential 
rotation might lead to a systematic error in the measurement 
of the projected obliquity. For a well-aligned system (A w 0°) 
the effect is much less important because the planet probes 
only a small range of stellar latitudes over the course of the 
transit (see also Gaudi & Winn 2007). 

We have limited empirical knowledge of differential rota- 
tion profiles in stars other than the Sun, and therefore we have 
limited ability to predict the impact of differential rotation on 
our observation of WASP-7. One relevant study is by Rein- 
ers & Schmitt (2003), who examined differential rotation in 
slowly-rotating stars (vsiniV < 20kms"'). In their sample, all 
stars having a logR'^j^ index between -4.80 and 4.65 showed 
signs of differential rotation, while more active stars had no 
discernible differential rotation. Their sample did not include 
any quieter stars {logR'j^f^ < -4.80). In this context, we should 
expect significant differential rotation for WASP-7, which ro- 
tates more slowly then 20 km s"' (Hellier et al. 2009), and has 
a low activity index logR'^j. of -4.981 (H. Knutson, priv. com- 
munication, 2011). We adopt the parameterization used by 
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Fig. 4. — Results for vsin/* and A, based on our MCMC analysis in the 
WASP-7 system. The gray scale indicates the posterior probability den- 
sity, marginalized over all other parameters. The contours represent the 2-D 
68.3%, 95%, and 99.73% confidence limits. The one-dimensional marginal- 
ized distributions are shown on the sides of the contour plot. 

Reiners Sl Schmitt (2003) and others, 

O(0 = f^eq(l-asin^/). (1) 

Here / denotes stellar latitude and ilgq the angular rotation 
speed at the equator, a is a dimensionless parameter which 
describes the degree of differential rotation. An a of would 
indicates uniform rotation and an a of 0.2 corresponds to Sun- 
like differential rotation. See also Hirano et al. (2011a), who 
used the RM effect to set an upper limit on the stellar differ- 
ential rotation in the XO-3 system. 

3.5. Modeling the anomalous RV 

In order to develop a model for the anomalous RV produced 
by the RM effect, we must take into account two other ef- 
fects: (1) the stellar absorption lines recorded by the spectro- 
graph are further broadened by the PSF of the spectrograph; 
and (2) the measured RVs represent the output of a Doppler- 
measuring code, which is akin to finding the peak of a cross- 
correlation between a stellar template spectrum obtained out- 
side the transit, and a spectrum observed during the transit. 
Some researchers neglect the second point, and model the 
anomalous RV as the shift in the first moment of the model 
absorption line rather than as the shift in the peak of the cross- 
correlation function. Depending on the system parameters 
this can lead to systematic errors in the derived parameters 
(Hirano et al. 2010, 2011b). 

Hirano et al. (2011b) developed an analytical description 
for the shift in the cross-correlation peak as a function of the 
transit parameters, the stellar rotation velocity and obliquity, 
the microturbulent and macroturbulent velocities, the differ- 
ential rotation profile, and the PSF width of the spectrograph. 
We use their approach to model the obtained RVs during tran- 
sit and also include a model as described above for the con- 
vective blueshift. See Figure 2 for an illustration of how the 
above described effects change the line shape and the time 
variation of the RM effect. 

3.6. Other RV variation sources 

To this point we have only discussed changes in the RVs due 
to the transit. To successfully model the RM effect we also 
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TABLE 2 

Parameters of the WASP-7 system 



Parameter 


Values 


Parameters mainly derived from photometry 


Midtransit time [BJDtdb-2400000] 


55446.635 ± 0.0003 


Period. P [days] 


4.9546416 ± 0.0000035 


cos Jo 


0.05 ± 0.02 


Fractional stellar radius, Rt,/a 


109 4- 0.011 


Fractional planetary radius, ftp/R* 


0.096 ± 0.001 


U1+U2 


0.34 ± 0.08 


Parameters mainly derived from RVs 


Velocity offset, [m s~' ] 


34 ± 5 


Linear slope during transit night, a[ [ms"'] 


-95 ± 41 


Quadratic slope during transit night, 02 [ms^'/day] 


500 ± 310 


sin r't sin A [km s" ' ] 


3.7 ± 0.3 


\/ V sin ii, cos A [km s"' ] 


0.24 ± 0.56 


Macro turbulence parameter, ^ [kms"'] 


6.4 ± 1. 




0.92 ± 0.06 


differential rotation parameter, a 


0.45 ±0.11 


cos i* 


0.18 ± 0.43 



Indirectly derived parameters 



Orbital inclination, io [°] 


87.2 


± 


0.9 
1.2 


Full duration, T14 [hr] 


4.12 


± 


0.09 
0.06 


Ingress or egress duration, T12 [min] 


27 


± 


6 
9 


Projected stellar rotation speed, vsini, [kms"'] 


14 


± 


2 


Projected spin-orbit angle, A [° ] 


86 
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need to model the change in radial velocity due to the orbital 
motion of the star. In addition, as we will see in section 4, 
WASP-7 shows a high level of RV noise on a timescale of 
days, which will introduce trends in the RV over the course 
of the transit night. Therefore our model allows for RV trends 
that are linear or quadratic functions of time, in addition to the 
RM effect and orbital motion. The physical interpretation of 
these trends is discussed in section 4. 

3.7. Quantitative Analysis 

Now that all the ingredients of our model have been intro- 
duced, we describe the various parameters in detail. The Ke- 
plerian orbital motion of the star is specified by the period 
(P), the time of inferior conjunction (T^), the semi-amplitude 
of the projected stellar reflex motion (K^), and a velocity off- 
set (7). Initially we assume the orbit to be circular, as Hel- 
lier et al. (2009) found no sign of an orbital eccentricity (al- 
though see Section 4). We also allow for additional linear and 
quadratic RV trends (oi and 02) on the transit night, 

RV(t) = j+RVoMt(t)+RVRM(t)+ai*(t-to)+a2*(t-tQf . (2) 

Here RVorbit(f) and RVRM(f) represent the radial velocities 
caused by the orbital motion and the RM effect, to is a point 
in time near the middle of our observation sequence. Since 
allowing for these trends we have relinquished any power to 
constrain K^, with the data, so we fix at the value reported 
by HelUer et al. (2009). 

The amount of light blocked at any given phase of the tran- 
sit depends on the location of the chord of the planet's path 
over the stellar disk, which is parameterized by the cosine of 



the orbital inclination (cos/q), the radius of the star in units 
of the orbital semimajor axis {Ri,/a), and the radius of the 
planet relative to the stellar radius (Rp/R^,). We use quadratic 
limb darkening parameters and mi™ to parameterize the 
relative intensity of the stellar surface in the wavelength re- 
gion of 5000-6200 A, the spectral region in which the RVs are 
measured. We chose limb-darkening coefficients uirm = 0.3, 
uirm = 0.35, based on the tables of Claret (2004). We allowed 
u\rm + U2rm to Vary freely and held fixed u\rm~U2rm at the tabu- 
lated value of -0.05, since the difference is only weakly con- 
strained by the data (and in turn has little effect on the other 
parameters). 

We parametrize the projected spin-orbit angle and the pro- 
jected stellar rotation by the quantities v^TsmT^cosA and 
V^vsinZtSin A, rather than vsin/* and A. We do this because 
our chosen parameters are less strongly correlated (e.g. Al- 
brecht et al. 2011). Our macro turbulence model was that of 
Gray (2005), assuming equal surface fractions of radial and 
tangential velocities, with a macroturbulence parameter ( sub- 
ject to a Gaussian prior of 6.4 ± 1.0 kms"'. This represents 
our expectation for a star of WASP-7's effective temperature 
(Gray 1984). To model the convective blueshift we assumed 
an outward blueshift of 1 kms"' at all positions on the star 
(Shporer & Brown 2011, and references therein). We further 
include a, the parameter which governs the strength of differ- 
ential rotation, and the stellar inclination towards the observer 
(/*) as free parameters. As Reiners & Schmitt (2003) found 
no star with a-\/sin/* greater « 0.4 (see their Figure 16) we 
also impose a prior on a which is flat until 0.4 and then falls 
off as a Gaussian function with a standard deviation of 0.1. 

Finally, we must specify the width of a Gaussian function 
representing the width of the lines due to both microturbu- 
lence and the PSF of the spectrograph. We chose cr = 3 km s"' 
for this purpose. Also, to represent the natural broadening 
of the lines we used a Lorentzian function with a width of 
1 kms"' (see also Hirano et al. 2011b). 

Additional information on the transit geometry comes from 
a transit light curve recently obtained by Southworth et al. 
(2011). They made their de-trended light curve available via 
VizieR. They also gave an updated transit ephemeris for the 
system, which is derived from the WASP discovery data in 
combination with the new light curve. We used those updated 
results for P and as priors (see Equation 3). We also fitted 
the light curve simultaneously with the RVs in order to pin 
down cos/o, Ri,/a, and Rp/R^,. We used the algorithm from 
Mandel & Agol (2002) to model the light curve which was 
obtained in the Gunn / filter. From Claret (2004) we obtained 
Ml = 0.17 and M2 = 0.36. We allowed U1+U2 to vary freely, and 
held Ml -M2 fixed at the tabulated value of -0.19. 

To derive confidence intervals for the parameters we used 
the Markov Chain Monte Carlo algorithm. The likelihood was 
taken to be exp(x^/2), where was defined as: 
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Fig. 5. — Orbital solution with zero eccentricity for WASP-7. The upper 
panel shows RV observations as function of the orbital phase, with periastron 
at zero phase. The CORALIE RVs are indicated by filled symbols and the 
HARPS data are shown by open symbols. The lower panel shows the residu- 
als between data and best-fitting model. 



The first two terms are sums-of-squares over the residuals 
between the observed (o) and calculated (c) values of the ra- 
dial velocity (RV) and relative flux (F). The following terms 
represent priors on some parameters, as mentioned above. Be- 
fore starting the chain we also added 10 ms"' in quadrature to 
the uncertainty of the PFS RVs to obtain a reduced ^ close to 
unity. In making this step we assumed that the uncertainties 
in the RV measurements are uncorrected and Gaussian. 

Our results are presented in Table 2. The best fits to the RVs 
and photometry are shown in Figures 1 and 3. Figure 4 shows 
the 2-d posterior density distribution for the two parameters 
of greatest interest for our study: A and vsin/*. 

For the parameters governed mainly by the photometric 
data, our results are consistent with those obtained by South- 
worth et al. (2011). We will discuss the reliability of some 
of the parameters found by the RV data, in particular spectral 
limb darkening, the evidence for differential rotation, and the 
stellar inclination in section 5.1, after investigating possible 
reasons for the excess noise in the RV data in the following 
section. 

4. STELLAR JITTER AND ORBITAL ECCENTRICITY 

The out-of-transit RV gradient measured with the PFS is 
steeper than would be expected from the previously published 
spectroscopic orbit. That slope, and the known orbital period, 
imply an orbital velocity semiamplitude of « 200 m s"' , in 
strong contrast to the published value of 97 ± 13 ms"'. (See 
the upper panel in Figure 1.) 

Archival RV data from the CORALIE and HARPS spectro- 
graphs, obtained during various orbital phases, show a large 
scatter around the best orbital solution in excess of the mea- 
surement uncertainties. This type of excess noise is com- 
monly referred to as "stellar jitter". In the following we will 
shortly investigate the timescale over which the RV scatter is 
correlated and test the possibility that the orbit is actually ec- 
centric. 

4.1. Excess noise 

We turn first to the previously reported RVs. Eleven mea- 
surements were obtained with CORALIE by Hellier et al. 



10- 
8 

6- 
4- 
2 

10- 




JD.QOl 



JDDl 
0.05 




10 

period [days] 



100 



Fig. 6. — Lomb-Scargle periodograms of the out-of-transit RVs. The 

upper panel shows the power at different periods in the RV data of the two 
data sets from CORALIE and HARPS. The highest peak occurs at the orbital 
period. The lower panel shows the power in the RV data after the best fitting 
model was subtracted. Significant peaks occur at a period of one day, and 
its higher harmonics. Most likely this is a consequence of the diurnal time 
sampling of the RV observations. 



(2009) and another eleven RVs were measured with HARPS 
by Pont et al. (201 1). We fitted an orbital model to both data 
sets separately to determine the root mean square (rms) resid- 
ual between the data and the best-fitting model. A high rms 
would indicate that our model neglects an effect which influ- 
ences the RVs, which for example could be activity of the 
star itself. We find that the scatter of both the CORALIE 
and HARPS RVs is not only greater than what these instru- 
ments normally achieve on bright stars, also the rms of the 
HARPS data (33 ms"') is slightly higher than the rms of the 
CORALIE RVs (31 ms"'). This is noteworthy as HARPS 
should achieve a greater precision then CORALIE, if photon 
noise is the limiting factor, as HARPS operates in conjunction 
with a 3.6 m telescope and CORALIE with a 1 .2 m telescope. 
This suggests that an additional source of RV variations may 
be present, which dominates the noise budget. 

Next we fitted each dataset individually, assuming a circular 
orbit and adding a term in quadrature to the internal uncertain- 
ties to produce a reduced of unity. In making this step we 
assumed that the errors of the RV measurements are uncor- 
rected and Gaussian. We needed to add 25 m s"' in quadra- 
ture to the internal uncertainties of the eleven CORALIE dat- 
apoints and 33 m s"' to the uncertainties of the eleven HARPS 
RVs. 

Using those inflated uncertainties, we fitted both data sets 
together, also assuming a circular orbit. Our fitting statistic 
was 



RV,(o)-RV,(c) 

0"RV,i 

/'-4;'9546416y 
0:'0000035 / 



rc.BJD- 2455446.63493 Y 
0.000030 / 



(4) 



making use of the ephemeris from Southworth et al. (2011). 

Figure 5 shows the phase-folded RV data. As expected 
from the previous fits, the rms residual is 31.4 ms"'. To 
understand this noise source it would be important to learn 
the timescale over which the RV noise is correlated. Looking 
at the Lomb-Scargle periodogram of the RV data, we find the 
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Fig. 7. — Constraints on the orbital eccentricity (e) and orientation (uj) 
of WASP-7 from HARPS and CORALIE data. The gray scale indicates 
the posterior probability density, marginalized over all other parameters. The 
contours represent the 2-d 68.3%, 95%. and 99.73% confidence limits. The 
eccentricity is constrained to low values. For an orientation of tj |90|° 
somewhat higher values of the eccentricity are allowed. That is particularly 
true for positive value of uj. The 1-d marginalized distributions are shown on 
the sides of the contour plot. 

most dominant peak at the orbital period of WASP-7b (Fig- 
ure 6, upper panel). After subtracting the best-fitting circular- 
orbit model, the periodogram of the residuals shows a strong 
peak at a period of one day and its harmonics (Figure 6, lower 
panel). This is most likely a result of the timing of the ob- 
servations. Many observations were obtained on consecutive 
nights. Apart from this pattern, the periodogram is not very 
informative. An autocorrelation plot does not reveal addi- 
tional information. 

The timing is different for the PFS observations. We ob- 
tained 37 data points during an interval of 7.5 hr The rms 
around our solution including the orbital model, the linear 
and quadratic acceleration, and the RM effect is 1 1 m s"' (Fig- 
ure 1). This relatively low scatter indicates that the correlation 
time of the RV variations is longer then a few hours. 

We note that WASP-7 is not the only early-type planet host 
star for which a large stellar jitter has been observed. Recently 
Hartman et al. (2011) found that the two planet host stars in 
the HAT-P-32 and HAT-P-33 systems have stellar jitter of « 
80 and 55 ms"'. 

4.2. Orbital eccentricity 

One possible contributing factor to the excess RV noise is 
that the orbit is actually eccentric, in violation of our mod- 
eling assumption of a circular orbit. To investigate this pos- 
sibility we repeated the MCMC analysis but allowed the or- 
bital eccentricity (e) and the argument of periastron (oj) to be 
free parameters. Our stepping parameters were ^/esinw and 
y/ecosLj, which have less correlated uncertainties for small 
eccentricity, and which correspond to a flat prior in e. 

Using the likelihood based on Equation 4, we derived the 
posterior probability, which is displayed in the uj-e plane in 
Figure 7. The rms for the best fitting eccentric solution is with 
30.9 m s"', only moderately smaller than for the circular-orbit 
model. We can see from Figure 7 that there is no clear de- 
tection of an eccentric orbit and that solutions with e > 0.2 
are generally disfavored except for values of oj near 90°, 
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Fig. 8. — Results for orbital eccentricity using synthetic data. Similar 
to Figure 7 but this time for an MCMC analysis of a synthetic data set which 
was derived from a circular-orbit model. The results indicate that higher ec- 
centricities are allowed for cu near 90° or —90° . 



for which somewhat larger eccentricities are allowed. Such 
an argument of periastron would indicate for transiting sys- 
tems a semi-major axis closely aligned with the line of sight 
(LOS). Interestingly, such an orbital orientation would lead to 
a steeper RV slope during the transit night, which is indeed 
what was observed with the PFS. 

However, the peak in the probability density near uj « 90° 
should not be taken as evidence that the orbit really does have 
this orientation. This is probably just the result of the fact 
that RV studies are better at constraining y/e x cosw than 
y/e X sinoj, and consequently we expect the confidence in- 
terval for e to be larger for u close to |90|°. Such an orbital 
configuration would lead to symmetric RV curves even for 
an eccentric orbit, i.e., the RV amplitudes at the quadratures 
would not differ. 

To investigate this point further we conducted numerical 
experiments similar to those carried out by Laughlin et al. 
(2005) for the HD 209458 system. We created simulated RV 
datasets with the same time stamps as the CORALIE and 
HARPS data sets. For this we assumed a circular orbit with 
the parameters of the WASP-7 system. We then added Gaus- 
sian perturbations to the model RVs, adopting a Icr uncer- 
tainty of 33 ms"', and then used the same MCMC analysis 
for these mock data as was used for the real data sets. A typ- 
ical 2-d posterior resulting from this experiment is shown in 
Figure 8. 

We found that greater values of e are permitted for orbital 
orientations of w |90|°. This should make us suspicious of 
any low-SNR detection of an orbital eccentricity with an uj 
of w |90|°. We conclude that the out-of-transit RVs give no 
indication of an eccentric orbit for WASP-7. This is in line 
with the upper limit e < 0.25 found by Pont et al. (201 1). We 
note that the w-dependent sensitivity to e is very similar to 
the V sin /^-dependent sensitivity to A that was explicated by 
Albrecht et al. (201 1), for low-SNR studies of the RM effect. 
In that case, higher values of the stellar rotation rate (v sin /) 
are allowed for A w 0° and w 180°, even when the underlying 
signal has no RM effect at all. 

5. DISCUSSION 
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5.1. Differential rotation parameters 

Our model took into account the possible effects of differ- 
ential rotation, through the parameters a and i^,. We found 
that the upper boundary of the posterior distribution for a 
is determined mainly by our prior. Specifically we found 
a = 0.45 ±0.11, with the prior enforcing a < 0.5 (see Ta- 
ble 2). If instead no prior is placed on a, then we found that 
the differential rotation parameter increases to values as high 
as 0.9, much higher than would be expected based on theory 
and on observations of other stars. 

This should make us suspicious. For |A| close to 90°, as 
is the case here, a parameter degeneracy exists between the 
limb-darkening profile and the degree of differential rotation, 
because both of those phenomena produce changes to the RM 
effect that are symmetric in time about the transit midpoint. 
Perhaps our assumptions are mistaken regarding the stellar 
limb darkening within the effective observing bandpass. For 
RM observations the effective observing bandpass is compli- 
cated to describe, depending as it does upon the density of I2 
absorption lines. Furthermore, there may be systematic dif- 
ferences between the limb darkening as observed in different 
absorption lines, and even within a single strong line, as dif- 
ferent parts of the absorption lines form in different depths 
of the stellar photosphere. Nevertheless, the fitted value of 
the center-to-limb variation (Mi„„ + M2/m = 0.92 ±0.06) is al- 
ready quite high; lowering this value would only cause a to 
converge to larger values. 

The data seem to be demanding a greater variation of the 
RM effect between the second and third contact than is deliv- 
ered by the expected amounts of limb darkening and differen- 
tial rotation. Some possible explanations are: 

• The stellar jitter exhibits correlations on the timescale 
of minutes to hours, which is not taken into account in 
our model. By fitting a linear and quadratic trend to the 
RVs, we have only taken into account correlations on 
longer timescales. 

• If WASP-7 is pulsating then velocity fields on the stellar 
surface, due to the pulsations, might lead to a different 
shape of the stellar absorption line than is created by our 
model, which does not include pulsations. The planet 
would occult during its transit different velocity com- 
ponents than expected. For this effect to be important 
the velocity fields do not need to change on a timescale 
comparable to the transit. 

• Another possibility is that our model of the convective 
blueshift is not correct. For hotter stars there exists a re- 
versed shape in the bisectors, although WASP-7 seems 
to be securely located on the 'cool' side of this 'granu- 
lation boundary' (see Fig. 17.18 in Gray 2005). 

With the present data we cannot determine if one of these 
or another effect is responsible for the relatively strong vari- 
ations in the shape of the RM signal. One should therefore 
view our particular results for the differential rotation, stellar 
inclination, and the spectroscopic limb darkening coefficients 
with some skepticism. Including these effects, or some ef- 
fects with a similar functional form, is however important for 
a realistic uncertainty estimation of vsin/^ and A as we will 
discuss in the next section. 

5.2. Stellar obliquity 
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Fig. 9. — Dependency of A on r*. The gray scale indicates tlie poste- 
rior probability density, marginalized over all other parameters. The con- 
tours represent the 2-d 68.3%, 95%, and 99.73% confidence limits. The 1-d 
marginalized distributions are shown on the sides of the contour plot. One 
can see the strong correlation between A and i*. This congelation vanishes if 
no differential rotation is present. 

As expected from the qualitative discussion we find A = 
86 ±6°, which is consistent with 90°. The projected stellar 
spin axis is lying nearly within the orbital plane. 

We did not use the measurement of the projected rotational 
velocity (vsin;,t = 17±2 kms"') by Hellier et al. (2009) as a 
prior constraint. This is because it is not clear what value of 
macroturbulence was assumed by those authors, and also be- 
cause if differential rotation is present then a systematic error 
could be introduced. 

The lack of knowledge on differential rotation and the stel- 
lar inclination did lead to an increased confidence interval for 
A. This can be seen from the Figures 4 and 9. When the stel- 
lar inclination departs from 90°, the planet covers higher stel- 
lar latitudes, which have decreased rotational velocities (and 
therefore a decreased RM effect) either at the beginning or 
end of the transit. This forces a higher vsin/^ and a change 
in A to compensate for the asymmetry in the transit RV curve. 
This degeneracy could be broken if we would have indepen- 
dent information on the stellar inclination. 

To estimate sin/^, we can use the technique of Schlauf- 
man (2010), involving a comparison of the measured value of 
V sin ii, with the expected value of v for a star of the given age 
and mass. Schlaufman found for WASP-7 a rotation statistic 
6 = -3.4, indicating that WASP-7 rotates faster then expected 
for its age and mass. A 8 near would indicate that the mea- 
sured vsin;\t is consistent with the expected rotation speed v 
for a star of a given mass and age. A G larger then would 
indicate an inclination of the stellar spin axis towards the ob- 
server 

We repeat his analysis with the new values for vsiniV, found 
here using the RM effect (Table 2), and the new mass, radius 
and age values from Southworth et al. (2011). With these we 
obtain an expected v of 11 ±4 kms"' and a rotation statistic 
6 of -0.8. This indicates that the projected rotation speed 
is consistent with a sini^, of 1 There is no indication of an 
inclination of the stellar spin axis towards the observer. 

This knowledge could be used as prior knowledge, or more 
correctly used in the model itself to decrease the uncertainty 
in the projected obliquity. However as WASP-7 is at the up- 
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per end of the mass range for which Schlaufman (2010) cal- 
culated his rotation, mass, age relationship and because of the 
complications due to differential rotation for which his rela- 
tion was not calibrated, we decided not to use this knowledge 
to reduce the uncertainty in A. 

5.3. Stellar jitter 

Analyzing the bisectors of the obtained spectra might lead 
to a reduced scatter, as was done for HAT-P-33 by Hartman 
et al. (2011). In particular the data taken during the transit 
night might be informative as the change in the bisectors will 
be correlated over the course of the night. 

If during some nights several RVs would have been ob- 
tained then the uncertainty in the orbital parameters of WASP- 
7 could be reduced in a similar approach to that used by 
Hatzes et al. (2010) for CoRoT-7. They used data obtained 
during one night to constrain a part of the orbit and allowed 
for a drift between different nights. However this approach re- 
quires a substantial amount of data and its success depends on 
the timescale over which the jitter is correlated, with shorter 
correlation timescales being advantageous. 

As we are mainly interested in the systems obliquity we 
only note here that the relation by Saar et al. (2003) employ- 
ing a correlation between vsin;^ and stellar jitter, leads to an 
expected jitter of 34 ms"', similar to the measured value. 

6. SUMMARY OF CONCLUSIONS 

We find that in the WASP-7 system the stellar spin axis is 
strongly misaligned with the planet's orbital axis, by 86 ±6° 
as projected on the sky. This observation strengthens the cor- 
relation found by Winn et al. (2010) and lends support to the 
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idea that systems with close giant planets generally started out 
with a very broad range of obliquities, and that the observed 
low obliquities of many systems are a consequence of tidal 
dissipation. 

Differential rotation and its imprint on the RM effect holds 
the promise of measuring not only the projections of stellar 
obliquities, but also the stellar inclinations. However, with 
the current measurement precision and uncertainties in other 
parameters such as limb darkening, no secure detection of 
differential rotation in WASP-7 can be made. We originally 
thought that WASP-7 might present a good testbed to search 
for differential rotation via the RM effect, as the misalignment 
of 90° maximises the RM signal originating from differential 
rotation. However for this angle the signal from differential 
rotation is also strongly correlated with stellar limb darken- 
ing. In addition WASP-7 displays a high degree of stellar jit- 
ter Therefore a system with a quiet star and a more moderate 
misalignment might be better suited to search for signs of dif- 
ferential rotation in the RM signal. 
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